Distributional Matter Tensors in Relativity 
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This paper uses products of distributions to obtain new junction conditions for relativistic shocks. 
In general, the shock is accompanied by a surface layer, and the new conditions generalize both 
Taub's jump conditions for shocks, and those of Israel and Kuchar for surface layers. In the non- 
relativistic limit, the surface layer is present only when the fluid is viscous or thermally conducting — 
a situation where the classical Rankine-Hugoniot conditions do not apply. Thus, our conditions 
properly extend all previous conditions, and provide complete Cauchy data needed to solve the full 
Navier-Stokes equations downstream of the shock. Since the associative law fails for our product, 
the residual uncertainty regarding the association of factors must be eliminated empirically. This is 
equivalent to fixing the correct initial form of the equations (such as the "conservation form" in the 
Euler case). 



I. INTRODUCTION 

Distributional matter tensors arise in problems con- 
cerning shells of matter, shock waves, line sources, and 
possibly some 'singularities' Such matter tensors 

involve problems with products of distributions. Here we 
give a quick overview of the theory of such products, 
and point out conditions under which such matter ten- 
sors may be regarded as distributional limits of smooth 
matter tensors. 

Distributional matter tensors are believed to be no 
more than simplified models representing a situation 
which is actually smooth; however, this belief is not al- 
ways justified. For example, for perfect fluids, one may 
arrange initial conditions so that (classically) the charac- 
teristics intersect, and genuine discontinuities develop 
strings are believed to be represented by 2-dimensional 
surfaces, and 'singularities' might also indicate a break- 
down of smoothness. 

The geometrical theory of distributions on a manifold 
is straightforward [5| except for the analytical problem of 
products. According to Taub 0], "Fortunately the prod- 
uct of such distributions (as arise) is quite tractable." 
Thus, if 9 is the Heaviside function, 9 2 = would sug- 
gest 29 ■ 9' — 9' , i.e., 9 ■ 8 = ^8. This is simple enough, 
except that 9 3 — 9 would suggest 39 2 ■ 9' = 9' which 
should be the same as 9 ■ 8 = |#. More generally, since 
x^ 1 (xS) = 0^8 = (x~ 1 x)8, the Schwartz impossibility 
theorem asserts that there does not exist an associative 
differential algebra ADD' such that the product in A 
agrees with the Schwartz product or the product of C°° 



"This paper was presented at MG5 (1988) and originally published 
in Proceedings of the Fifth Marcel Grossmann Meeting on General 
Relativity, ed. D. G. Blair and M. J. Buckingham, World Scientific, 
Singapore, 1989, pp 419—22. Apart from the practical value of the 
new (and still unused) shock conditions, the paper illustrates an 
important point in the new philosophy of mathematics proposed 
in my recent book, Cultural Foundations of Mathematics, Pearson 
Longman, 2007. Hence, it has been reproduced verbatim on the 
archive (with the addition of two footnotes and the abstract). 



functions. 

In contrast to Taub's naive approach, Parker's sug- 
gestion to use Hormander's product would exclude even 
the theory of shells, since 9 ■ 8 is not defined then. At the 
other extreme, Geroch and Traschen have proposed a no- 
tion of a regular metric (including possibly discontinuous 
metrics) in such a way as to exclude entirely products of 
distributions. This seems too strict, since such products 
have been freely used in renormalization theory and are 
known to 'work' ! Such regularity conditions also depend 
on the precise form of the equations used, since p.d.e.'s 
with equivalent smooth solutions may have incquivalent 
weak solutions. 



II. PRODUCTS OF DISTRIBUTIONS 

Currently, the real problem is to select one definition 
from amongst the many that are available. According to 
the Fourier transform method @, 0, H, Q , the product is 
essentially defined for /, g £ D' by / -g = (/ ® g), when- 
ever the r.h.s. exists, ® being convolution. This prod- 
uct is Tol Ilia particular case of the sequential product 
[H [rMJUl f-g = D'- lim(/ ® 8 n ) ■ (g ® p m ) when- 
ever the r.h.s. exists for 'all' <5-sequences 8 n , p m . Putting 
8 n = Pm = an approximate identity, and requiring the 



limit to be independent of the choice, we get the sequen- 
tial 'model' product 16]. Tighter conditions on 8 n would 
allow more products to exist, but could lead to a sit- 
uation where physically acceptable regularisations have 
non-unique limits. 

Colombeau's [l?], HH associative differential algebra 
GD D', and the product does not agree with the product 
of C°° functions. The product of any f,g£ D' always 
exists in G but admits an 'associated' distribution iff [19[ 
the sequential model product exists. However, the prod- 
uct is not coherent with association so that 8 ■ 8 need 
not have a unique associated distribution. This prod- 
uct is thus analogous [20] to the product in the non- 
standard space *E =* C°°. A similar ambiguity exists 
for products defined by the Hahn-Banach method used 
in renormalization theory. There one defines [2l|, [22j 
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d a (f ® g) = /(g) d a g, provided the r.h.s. exists. The 
Hahn-Banach extension to the whole space is non-unique 
upto arbitrary linear combinations of 6 and its derivatives 
to order a. Thus 6-5 = AS, 6 2 = B8, with A, B arbitrary 
constants. 



I define [23| / • g = \ 
* denotes the non-standard extension to 



{(f®8 u )g 



(g ® S u )f}, where 
D' (i.e., *E) 

and a; is a positive infinite integer. The Leibniz rule 
holds, but the associative law fails. The product of any 
/, g £ D' always exists, is unique and coincides with 
the sequential product when the latter exists. We have 
9 ■ 5 = i<5, and 5 2 = 5^(0)5, where <5 W (0) is infinite. Nev- 
ertheless, for applications to weak solutions of p.d.e.'s 
a kind of linear independence [23, [24( ensures that the 
final results are standard. For applications to renormal- 
ization, it may be shown that the infinities arising in this 
way actually do cancel [25] . 



M - o, 

p=(2 V + X)[v], e=~[T\, 



(pv)\[w + -v 2 



[p + pv 2 ] = (2tj + A) [«,„], 
(2r) + \)[vv, n ] +k[T,t 



(3.1') 
(3.2') 

(3.3') 

e,t, (3-4') 



where v is the fluid velocity in the direction of the normal 
n to the shock, terms with a ~ are defined as before, w is 
the enthalpy, and the remaining notations are standard. 
Thus (|3.2[) is non-trivial. 

(|3. 1[) — (|3.4[) are still insufficient to determine adequate 
Cauchy data on E, since an arbitrary equation of state 
may not be used on E, but, using the first law of ther- 
modynamics, and S 2 , one gets 



III. SHOCK WAVES IN REAL CONTINUA 

For reasons of space and time, I will mention here only 
applications to shock waves in real continua (e.g. in fluids 
with viscosity). The problem here is not that of using a 
smooth solution of the viscous-fluid equations to interpo- 
late a shock solution in a perfect fluid (as e.g. in Ref. 4). 
Rather the problem is to develop solutions of the viscous 
equations into the region downstream of the shock. For 
this purpose the classical conditions of Taub or Rankinc- 
Hugoniot fail to provide adequate Cauchy data on the 
shock hypersurface E, and the standard interpolation 
techniques are also useless. 

For timelike E, the jump conditions across E are 



[Nu 1 ] = 0, 
f Ql = 0, 



-[T% 
-IT 11 ]- 



(3.1) 
(3.2) 
(3.3) 
(3.4) 



Here T = lim < T, cj) > as supp (f> | E, for <j) 6 D 
with <f> = 1 on a neighbourhood of E, T being the stress 
tensor of the continuum, u denotes covariant differenti- 
ation w.r.t. the intrinsic geometry of E, [] denotes the 
jump, and | the mean value. (I am not convinced that 
the stress tensor can be defined without using admissible 
coordinates: x 1 is the gaussian coordinate normal to E.) 

Note that a surface layer is present, and (|3.2p - (|3.4|) 
are the conditions for surface layers except that (13.211 is 
identically satisfied in the formalism of Israel [2fl [27| 
and Kuchar (28|. For a perfect fluid, all terms with a 
vanish and the conditions reduce to the Taub conditions. 
But for a viscous and thermally conducting fluid, the 
nonrelativistic limit is 



\l4 = jfi\p] + T\h 1 [S\, 



p : 



(NT)[N] _ 
[NT] fI 



ho 



(3.5) 
(3.6) 



where ho and hi are uncertainty factors due to the failure 
of the associative law, 1 S is the entropy, and p is the 
chemical potential. ()3. ![) — (|3.6[) suffice to fix Cauchy data 
on E. Incidentally, the weak-shock non-relativistic limit 
of KB is 



T\[S] = 



W\j 



■[pf + {2r 1 + \)v, n \[v] 



K 1 

-[T,n] ~ -e,t 
J J 



(3.5') 

giving an exact expression for the entropy change. 

The distributional matter tensors may also be thought 
of as simplified models: the weak solutions g a p € C° are 
the uniform limit on compacta of the smooth g a p ® S n 
which satisfy Einstein's equations with the matter tensor 

MOn = 87tT ^ ® <U - g a- ® <*»• (3-7) 



(T;j n -> T M „ in D' since G'^ -> in (E2D, provided 
only that the products of distributions appearing in Ein- 
stein's equations exist according to the sequential model 
product. A background of positive energy is a sufficient 
condition for positivity of energy with this T'. If asso- 
ciativity is taken for granted, and products such as 9-8-d 



1 This is a desirable feature of the theory, since the failure of the 
associative law explains why p.d.e.'s with equivalent smooth solu- 
tions may have inequivalent discontinuous solutions — a counter- 
intuitive point on which Riemann slipped (he incorrectly sup- 
posed that energy and entropy conservation were equivalent 
across a shock). This simple explanation is not available with 
associative products like Colombeaus's, and an earlier attempt 
(Raju and Colombeau 1987, unpublished) to compare this prod- 
uct with Colombeau's by checking possible differences in the 
resulting junction conditions (for the non-relativistic case) was 
abandoned. 
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arise, different regularisations would have different lim- 
its. On the other hand, it may be possible to make sense 
out of a wider class of genuine line sources if my product 
is used, and a definite association of factors can be fixed. 
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